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Abstract 

Popov  type  frequency  domain  conditions  for 
stability  of  feedback  systems  containing  ferromagnetic 
hysteresis  non-linearity  are  established. 

I.  Introduction 

Popov  criterion  and  its  extensions  consider 
non-linear  elements  that  are  memoryless  and  pass 
through  the  origin,  i.  e. ,  f(0)  *  0,  For  an  important 
class  of  non-linearities,  ferromagnetic  hysteresis, 
none  of  the  above  conditions  are  satisfied,  i.  e. ,  it 
is  neither  non-dynamic  nor  pass  through  the  origin. 
Therefore,  to  analyze  the  stability  of  systems  con¬ 
taining  this  type  of  non-tinearity.  appropriate 
modifications  to  Popov's  approach  should  be  made. 

Published  material  to  tackle  this  problem  is 
scarce.  The  only  work  known  to  us  is  by  Lecoq  and 
Hopkin  fl],  where  by  letting  the  derivative  of  their 
input  signals  to  belong  to  exponentially  weighted  Lz 

spaces,  they  obtained  bounded  input- bounded  output 
stability  for  systems  containing  hysteresis  non- 
linearities. 

In  the  present  paper  we  analyze  the  stability  of 
feedback  systems  of  the  form  shown  in  fij.  (la), 
where  N  is  a  ferromagnetic  hysteresis  non-linearity 
and  H  is  a  linear  element.  The  analysis  is  done  by 
substitution  of  the  model  for  the  hysteresis  proposed 
by  Chua  and  Stromsmoe  [2].  Then  the  concept  of 
passivity  is  utilized  to  derive  Popov  type  frequency 
domain  conditions  on  the  linear  element  H  for  stability 
of  the  feedback  system.  It  will  be  shown  that  if  the 
same  conditions  as  in  the  classical  Popov  criterion1 
are  satisfied  by  inputs  Uj  and  U2  and  linear  element 

H  then  the  feedback  system  of  fig,  (la)  is  stable  if 
the  non-lineaT  element  N  is  a  ferromagnetic  hystere¬ 
sis. 

11.  Hysteresis  Modeling 

The  model  for  ferromagnetic  hysteresis  of  Chua 
and  Stromsmoe  [2]  is  given  by 

•  *o[x(t)  -  f  o(y(t))]  (2.  I)2 

where  x(t)  and  y(t)  are  real- valued,  continuous  input 
and  output  signals  of  the  hysteresis  non-linearity 

'Research  supported  in  part  by  AFOSR  Grant  80-0013 
and  in  part  by  NSF  Grant  1NT-8302754, 

1  By  the  classical  Popov  criterion,  we  refer  to  the 
earliest  version  of  the  result  obtained  by  Popov  j3] 
and  also  derived  by  different  approaches  ir.  14],  bi 
(6],  and  not  later  generalizations  of  the  result  by 
Yakubovich  [7]  among  others. 

^o:  Functional  Composition. 
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Table  1 
Notations 

Meaning 

Field  of  real  and  positive  real  numbers. 
The  space  of  signals  such  that  r  |x(:>  ,dt 
exists.  -  = 

The  space  of  bounded  signals. 

The  space  of  signals  which  are  square 
integrabte  on  every  bounded  interval 
f0,  T  ]  9], 

The  Hilbert  space  of  sipna.s  writer,  art- 
square  integrable  on  (-*  , -)  with  ir.'.er 
product  <  x,  y  >. 

r“  y*(t)  x(t)  dt 

•  m 

J  <  X,  X  > 

rV<t)  x(t)  at 

o 

/<  X,  x>_ 


Truncated  x  |b] 

Set  of  instances  of  time  of  interest. 
Convolution  Algebra  ]. 

Real  part  of  a  complex  quantity. 
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representing  th*  current  i(t)  and  the  flux  linkage  c(t) 
•f  an  inductor  (transformer);  and  g  and  f  are  strictly 
monotonies Uy  increasing,  differentiable,  onto  func- 
tions  enjoying  the  important  property  of 


t  +  t  >  0  (3.  4) 

Then  the  feedback  system  is  finite  gain  L^- stable. 
Proof:  See  for  example  [6],  g 


g(0)  •  f(0)  >  0 


IV.  Main  Results 


Equation  (2. 1)  models  the  behavior  of  ferro¬ 
magnetic  hysteresis  successfully  and  with  very  good 
accuracy.  It  predicts  the  expansion  of  the  area  of 
the  hysteresis  loop  with  increasing  frequency  and 
or  edicts  minor  hysteresis  loops  such  as  commonly 
•ccur  when  a  d- c  plus  periodic  input  is  applied. 


After  plotting  the  hysteresis  loop  for  a  conven¬ 
ient  signal,  simple  procedures  are  given  in  [2]  to 
determine  g  and  f  for  that  loop.  When  non-linear 
functions  g  and  f  are  determined,  they  can  be  substi¬ 
tuted  la  the  model  of  equation  (2. 1)  to  predict,  with 
good  accuracy,  the  hysteresis  shape  and/or  its 
output  for  any  arbitrary  input.  For  further  detail 
and  examples  see  f2j. 

HI.  Passivity  and  Stability 


Definition  (3, 1)  [61:  Let  H:  -  L2e<  Then 

H  is  passive  iff  there  exists  some  constant  S€R  such 
that 

<  Hx,  x  >,j.  >5  V  x  €  L^^ 

V  T€J 

B 

Definition  (3.2)  (fel:  Let  H:  L^(  -L^.  Then 

H  is  strictly  passive  iff  there  exists  i  >  0  and  some 
constant  9€R  such  that 

<Hx.  x>T  >  JVxt'  +6  **€L2e 

VT€J  H 

Definition  (3.  3)  fill:  The  feedback  system  of 
fig.  (la)  is  said  to  be  finite  gain  L,- stable  if 

»)  «r*2’yl,y2'L2  V  Vu2€L2 


b)  There  exists  constants  pj  and  p2  such  that 

I'cj  I1, 1'  «2):,  J'  yj  P,  )*y 2  !•  £  !'  Uj|‘  ♦  Pjl  uj!  ^ur  u2  -L2 


In  the  following  well-known  theorem,  the 
concept  of  passivity  is  used  to  establish  finite  gain 
L2>stability  of  feedback  system  shown  in  fig.  (la), 

where  N  and  H  are  considered  to  be  operators  in  the 
general  sense. 


Theorem  (3. 1):  Consider  the  feedback  system 
shown  in  fig.  (la) 


u,  -  He2 


tt2  +  N*1 


Substitution  of  the  model  given  by  equation 
(2.  1)  for  the  hysteresis  non-linearity  N  gives  the 
feedback  system  of  fig.  (2).  Note  that  although  the 
standard  magnetic  hysteresis  non-linearity,  as 
shown  in  fig.  (lb),  is  the  plot  of  flux  linkage  (::(t)) 
vs.  current  (i(t))  of  an  inductor  (transformer),  but 
from  circuit  analysis  point  of  view,  the  input  and 
output  of  the  model  replaced  for  N  as  shown  in 
fig.  (2)  are  current  through  and  voltage  across  the 
inductor  (transformer).  * 

Next,  the  main  stability  result  is  presented. 

Theorem  (4.  1):  Consider  the  feedback  sytem 
of  fig.  (2),  where  h(t)€  Lj(R+)4,  and  h(t)f 

Assume  that  for  any  Uj,u2f  L2>  there  are  solutions 
*l,e2,yry2^^2e"  ^  4  constant  q>  0  exists  such 
that  for  some  constant  6 

Re  {(1+qju.)  H(ji);  =  6>0  V*>0  (4.1) 

Then  V  Uj.Uj.  u2*  L2 

a)  (i)  e1.e,.e2.y1.y2.y2€L2. 

(ii)  there  exists  constants  Pj  and  #uch  that 

,ei  ei  •  e2  -  yr-  v2l.!  y2 

CCjt/uj'  +q  u1)  +  C2'u2' 
i.e.  ,  finite  gain  L2-stability. 

b)  Cj.  e2.yj,y2^  L^  are  continuous,  and  goto  zero 
as  t  -  * . 

B 

Hysteresis  N 


Fig.  2 


where  H,  N  :  L2#  •*  l-2e.  Assume  that  for  any 


U|,u2€L2  there  are  solutions  ij,«2(  L^. 

Suppose 

that  there  are  real  constanta  v,  1,  and  t 

such  that 

r««:T  <  v!vt 

(3.  1) 

<  Hx,  *  >i>  >  6  V  x 

(3.  2) 

<  x.  Nx  >T  >  t  |?Nx 

(3.  3) 

V«(  L2e,  VT«J 


Proof:  See  appendix.  B 

Upper  and  lower  bounds  of  non-linearities  g 
and  f  can  be  taken  into  consideration  to  obtain  less 
conservative  classes  of  linear  element  H(s).  As  a 
example,  an  upper  bound  on  g  is  exploited  in  the 


^The  voltage  across  an  inductor  (transformer)  is 
proportional  to  rate  of  change  of  the  flux  linkage, 
constant  of  proportionality  being  the  number  of 
turns. 

*h(t)  I#  the  impulse  response  of  H(s). 


Under  these  conditions  if 


-'.iVA.-w- 


following  Corollary. 

Corollary  (4.  II:  Consider  the  feedback 
•yetem  of  fig.  2,  where  g?  sector  10,  k),  h(t)€  Lj(S^), 

li(t)€«4.  Aname  that  for  any  Uj,u2€L2>  there  exieti 
•olutione  ej,e2.yj.y2€  1^^.  If  a  constant  q  >  0 
exists  euch  that  for  some  constant  t 

Re  (0<fqj*)H(j*)l  +  £  *  C  >  0  V  a>  >  0  (4.  2) 

then  V  Uj, u j, u2?  l-2  conclusions  of  Theorem  (4.  1) 
hold.  a 

Proof;  For  outline  of  the  proof,  see  appendix. 

B 

V.  Conclusion 

Popov  type  frequency  domain  conditions  for 
stability  of  feedback  systems  containing  ferro¬ 
magnetic  hysteresis  non-linearity  are  established. 

To  obtain  the  results,  model  of  Chua  and  Stromsmoe 
[2]  for  hysteresis  is  employed  and  the  concept  of 
passivity  is  utilised. 


To  simplify  the  proof  of  Theorem  (4.  1)  and 
Corollary  (4.1),  the  following  two  lemmas  will  be 
proved  first. 

Lemma  (A.  I);  Let  q  >  0.  g€  sector  (0,  *  ). 
Then  the  system  of  fig.  (A.  I)  is  passive  (6). 


I 

i  l 

!■*  q 

g(* )  f 

Fig.  (A.  1) 

Proof:  cy,x>j  *  <  g(;),  o  +  qr>T 

*  <  g(-).  ?>T  +  q<  g(T).  r>T 

(i):  <  g(r).  3  >_  >  0  because  g€  sector  (0,  •) 

T 

(ii):  q<  g(r).  r  >T  »  q  r  g(r)  ?  dt 
1  o 

_c(T) 

*  q  T  8(3)  d; 

0(o) 

Define  G(r)  =  ?  g(i)d!  .  where  C;  R  -R. 

4  o 

Clearly  C(r)  >  0  V  x€R.  Then 
q<(l:).:>T>-qC[f(o)]  V  T  >  0 

(l)  and  (ii)  implies  <  y.x  >_  >  -q  G[  c(o)]  V  T  >  0 
_  v  x;"L,  . 

Passivity  follows.  ‘e 

B 

Lemma  (A.  2);  Le:  q  >  0.  I€  sector  (0.«), 

^  *  0.  Then  the  system  of  fig.  (A.  2)  is  passive. 


■jf(-)  I - 

Fig.  (A.  2) 


Proof:  <  y,  x  *  <  o  +  q?,  c  > ^ 

«  <f(:).e>I4qc  5r(f(r,1.  c- >T 

(i)  <f(o).e>T  > -F[c(0)]  V  T  >  0 

where  F:  R  —  R,  by  lemma  (A.  1.  (ii)). 

<U>  q<£[f(OU>T  «  ,  r0T^[f(=)l|f  dt 

•q m|fi2dt>o 

because  y;  [f(s)]  >  0  and  {—  ]2  >  0 
(l)  and  (ii)  implies  that  <  y.x>_  >  -F["(0)]  V  T  >  0 
Passivity  follows. 


ai+qij 


y2  e2  2 

Fig.  (A.  3) 

Proof  of  Theorem  (4.  1):  (a):  iiy  ir.clusio-.  of  the 
multiplier  (1+qs),  q  >  0.  transform  the  feedback 
system  of  fig.  (2)  to  the  one  shown  in  fig.  (A.  3). 

(i) :  <y,.e,>x  =  <y,.e>T+ <  y,.f>T 

(is):  <V|.e>T  >  -qCle(0)l  by  lemma  (A.  1) 

(ib):  <  y  j .  R  >j  >  -qF[c(0)]  by  lemma  (A.  2) 

Then  (ia)  and  (ib)  implies  that  feed  forward  block  is 

passive. 

(ii)  (l*qs)H(s)  is  strictly  passive  if  inequality 
(4.  1)  is  satisfied.  Since  h(c)-:Lj(F._)  ar.-i 

li(t)?e(.  (l-qs)H(s)  has  finite  pair.. 

Therefore,  by  Theorem  (3.  1)  the  feedback 
system  of  fig.  (A.  3)  is  finite  gain  L^tub'.t 

V  Uj.Uj.Uj,?  Lj,  i.  e. .  8j.e2,yj,y2-  L2  and 
V « j  1 .  )'*2  .  y,  .  y  2  £  k  J  V  qul  *:2  «2  • 

From  Fig.  (A.  3).  y2(t)  =  m(t)  *■  y 2 ( t )  where 

m(t)  =  L'*{  m{t)c  Lj  and  y2(:)'L2.  there¬ 

fore,  y2(t),  y2(t)?  L2  [6,  Appendix  C],  Furth.erm.  ire, 
l,y2(t)!  Vm(t)  *  y2(t)'  L  <  '  m(t)  L  y2(t)  ^ 

|6,  App.  C]. 

But  V  m(t)  .  is  finite  (=  Co:. star.:  C). 

*"1 

Therefore 
-  - 

L*,|.’l;  Inverse  Laplace  Transform. 


.  -  .  .N 


,V  -• 


y2<*>  1L2sci*2<*>,L2 

<  CCji  uj+qi,  \*  C^2Wl2 

£C6i^-i’V,,'i,»,V+C9i!'“2,L2 

Cn  the  other  hand,  y2(t)  *  m(t)  *  y2(t).  Then 

<  c-p^lujl^t qi'u,!!^)*  c-s^u2'^ 
beceuee  f  m(t)  ,  ie  finite  too. 

h 

Similarly.  «|(t)  *  m(t)  *  Sj(t).  Therefore, 
similar  conclusions  for  *j(t)  follow  Immediately. 
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(b):  and  Si.tjf  Lj  implies  that 

y,  e€La,  are  continuous.^and  go  to  zero  as  t  -•(11). 

Since  the  model,  i.e. ,  equation  (2. 1),  is  a  continuous 
mapping  from  input  to  output  [2],  therefore,  tjf  L#, 

and  Cj(t)  -0  as  t  —  •  implies  that  the  same  properties 
hold  for  y2(t),  i.e.,  yj(t)<UB,  is  continuous,  and  go 
to  zero  as  t  - 


(10)  MacFarlane,  A.  G.  J.  ,  Engineering  Systems 
Analysis.  Reading,  Mass.,  Addison-V  esley. 
1964. 

[11]  Willems,  J.  C. .  The  Analysis  of  Feedback 
Systems.  Cambridge,  Mass. ,  MIT  Press, 
1971. 


Similar  conclusions  for  e,  are  immediate. 

2  BS 

Proof  Outline  of  Corollary  (4.  1):  Apply  a 
positive  feedback  of  gain  around  g.  To  compensate 
for  it,  apply  a  positive  feed  forward  with  gain  1  /k  to 
K(s).  Let  g  *  (g~  ■([)”*•  Then  g€  sector  (0,  *)  and 

g(0)  =  0.  Following  the  same  procedure  as 
Theorem  (4.  1).  Conclusions  are  immediate. 
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